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ABSTRACT 



We present a result of cross-correlating the Infrared Astronomical Satellite Faint 
Source Catalogue (IRAS FSC) with the spectroscopic catalogues of galaxies in the 
. Fourth Data Release of Sloan Digital Sky Survey (SDSS), the Final Data Release 

O ■ of 2dF Galaxy Redshift Survey (2dFGRS) and the Second Data Release of 6dF 

(— i | Galaxy Survey (6dFGS). We have identified 324 ultraluminous infrared galaxies 

O (ULIRGs) including 190 newly discovered ULIRGs, and 2 hyperluminous infrared 

I 1 galaxies (HLIRGs). Adding these new ULIRGs, we increase the number of known 

ULIRGs by about 30 per cent. The reliability of the cross-correlation is estimated 
using the likelihood ratio method. The incompleteness of our sample introduced by 
the identification procedure in this study is estimated to be about 5 per cent. Our 
sample covers the redshift range of z=0. 037— 0.517 with a median redshift of z=0.223, 
which is larger than that (z=0.184) of the sample of previously known ULIRGs. 
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1 INTRODUCTION 

The interest in luminous infrared galaxies, in particular ul- 
traluminous infrared galaxiefl (ULIRGs) and the hyper- 
luminous infrared galaxie^] (HLIRGs), has been growing 
since the launch of Inf rared Astronomical Satellite (IRAS; 
iNeugebauer et all 1 19841 ) in 1983. A great deal of effort has 
been made to understand the origin of the enormous in- 
frared luminosities of these populations and the time evo- 
lution of individu al ULI RGs, which are summarised in 
ISanders fe Mirabell i|l996l ) and lLonsdale. Farrah. fc Smithl 
(2006). It is generally accepted that dust heated by 
some combination of starburst and active galactic nuclei 
(AGN) activity is responsible for the IR luminosity (e.g. 
iFarrah et all [2003L However, it is not yet clear which is 
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1 We define ULIRGs as the galaxies whose infrared (1-1000 pm) 
luminosities are greater than 10 12 Lq. 

2 We define HLIRGs as the galaxies whose infrared luminosities 
are greater than 10 13 Lq . 



the dominant power source. There are sever al suggestions 
for the time evolution of individual ULIRGs. I Sanders et al.l 
(1988a) suggested that ULIRGs may be an initial , dust- 
shrouded stage of optical QSOs, while IFarrah et all (|200lh 
proposed that ULIRGs are not a simple transition stage 
from galaxy mergers to QSOs but evolve through diverse 
paths. On the other hand, several authors suggested that 
ULIRGs might evolve into moderately massive (L* ) field 
ellipticals (iKormendv fc Sanders] Il99a iGenzel et ail l200ll : 
iTacconi et al.ll2002h . 

Infrared space telescopes have dri ven much recent 
progress (see lGenzel fc C csarskv 2000 and lVerma et al ]|2005l 
for a review). IGenzel et al.l ( 19981 ) showed that several di- 
agnostic lines in m id-infrared spectra using Infrared Space 
Observatory (ISO; iKessler et all Il996l ) are very powerful 
in characterizing the power source of ULIRGs. They con- 
cluded that most ULIRGs (80%) are powered by star- 
bursts, but at least half of their sample require both a star- 
burst and an AGN. Similar results are found inlLutz et al 



j 1998 1 ) and iRigopoulou et alj (|l999l ). Moreover. [Tran et al 
i|200ll ) showed that ULIRGs with L ir <10 12 4 L Q are mostly 
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starburst dominated, while ULIRGs with L; r >10 
more likely to contain A GN. The Spitzer Space 



Tele- 



scope l|Werner et al.l 120041 ) and submillimetre/millimetre- 
wave cameras such as SCUBA are now extending samples 
of ULIRGs to higher redshifts (z > 1). Several au thors 
i|Houck et alJ 120051 : iLutz et al.l 120051 : lYan et al.l [20051 ) took 
spectra of optically faint and infrared luminous population 
using Spitzer, detecting broad spectral features such as PAH 
(Polycyclic Aromatic Hydrocarbon) emission and silicate 
absorption. They showed that the majority of these pop- 
ulation are at high redshift (z ~ 2) and have mid-infrared 
spectral shapes similar to local AGN-dominated ULIRGs. 

Sin ce th e publication of the IRAS Point Source Cat- 
alogue (|l988l : hereafte r PSC) and the IR AS Faint Sources 
Catalogue - Version 2 l|Moshir et al.lll992l . hereafter FSC92), 
many wide-area redshift survey follow-up campaigns have 
been conducted. In total, these have led to a heteroge- 
neous compilation of a few hundred ULIRGs from sources 
such as th e IRAS 1.2 Jy Re dshift Survey l)Strauss et al.l 
1 19901 . 1 19921: iFisher et all 1 19951) the IRAS 1 Jy Survey of 
ULIRGs (|Kimfc Sanders! 1 19981). the QDOT a ll-sky IRAS 
galaxy redshift survey l)Lawrence et al.l 1 19991). the P oint 
Source Catalogue red shift survey (ISaunders et al.ll200Cl ). the 
FIRST/FSC sample [[Stanford et al.ll2000l). and the Revised 
Bright Galaxy SampldFl ( Sanders et al.ll2003l). The smal l area 
surveys are summarised in [Sanders fc Mirabell (| 1996T ) and 
lLonsdale. Farrah. fc Smithl i|2006l ). The majority of ULIRGs 
found in all sky redshift surveys are from the nearby universe 
(z < 0.5) and have abundant multi-wavelength data. Such 
objects are useful prototypes for the study of high redshift 
infrared luminous galaxies. However, these redshift surveys 
are still far from being a complete spectroscopic survey of 
all IRAS sources. Therefore, many ULIRGs remain undis- 
covered in FSC92. 

More recent galaxy redshift surveys such a s the 
Sloan Digital Sky Survey (SDSS, lYork et al l |2000|). the 
2dF Galaxy Redshift Survey (2dFGRS. IColless et af 



2001) and the 6dF Galaxy Survey (6dFGS, iJones et al 
2004 ) have provided redshifts for much larger samples 
of galaxies. These data will clearly be very useful for 
finding new ULIRGs in FSC9 2, wh ich has been at- 
tempted by several studies |Cao et alJ \ 
iPasquali. Kauffmann. fc Heckmanl 120051). However, 



i|2005l ) and IPasquali. Kauffmann. fc Heckmanl (|2005l ) cross- 
correlated the IRAS sources with the SDSS spectroscopic 
sample of galaxies using a circular matching tolerance in- 
stead of using a positional error ellipse of the IRAS sources. 
ICao et all (|2006l ) performed a similar cross-correlation using 
the formal positional error ellipses of the IRAS sources, but 
they used only the Second Data Release of SDSS. 

In this paper we present a search for ULIRGs and 
HLIRGs from a cross-correlation of FSC92 with the most re- 
cent spectroscopic catalogues of galaxies in SDSS, 2dFGRS 
and 6dFGS, using positional matches in the IRAS error el- 
lipses. The data used for the cross-correlation and the iden- 
tification algorithm for ULIRGs and HLIRGs are presented 
in Section 2. Our main results are presented in Section 3, 
and are discussed in Section 4. We summarize and conclude 



3 The earlier vers i on of this survey is the B right Galaxy Survey 
llSoifer et al.lll986l. Il987l. Il989l: ISanders et aTlll995l) . 



in Section 5. Throughout this paper, we assume that Ho 
70 km s" 1 Mpc" 1 , D. Q = 0.3 and A = 0.7. 



2 IDENTIFICATION OF ULIRGS AND 
HLIRGS 

2.1 The Samples and Cross-Correlation 

We cross-correlated the IRAS sources in FSC92 with 
the spectroscopic s ample of galaxies in the Fourth Data 
Release of SDSS (|Adelman-McCarthv et all 120051 . here- 
after SDSSDR4), the Final Data Release of 2dFGRS 
|Colless et all l200ll . hereafter 2dFGRSDRF) an d the Sec- 
ond Data Release of 6dFGS l|Jones et al 1 120051 . hereafter 
6dFGSDR2). The FSC92 contains about 173,000 infrared 
sources with measured fluxes at 12 u.m, 25 u.m, 60 u.m 
and 100 urn. The SDSSDR4 and the 2dFGRSDRF have a 
similar median depth of £=0.11. The SDSSDR4 contains 
about 473,000 galaxies over 4,800 deg 2 , and the 2dFGRS- 
DRF contains about 246,000 galaxies over 2,000 deg 2 . The 
6dFGSDR2 which contains about 89,000 galaxies, covers a 
wider area (~13,600 deg 2 ) than the other surveys, but has 
a shallower median depth of z=0.05. 

The positional uncertainties of the IRAS sources are 
different for the in-scan direction and the cross-scan direc- 
tion (typically 5" for the in-scan direction and 16" for the 
cross-scan direction, 1 er), and vary from source to source 
(1" — 13" for the in-scan direction, or minor axis of positional 
uncertainty ellipse and 3"— 55" for the cross-scan direction, 
or major axis of positional uncertainty ellipse, 1 a). The 
positional uncertainties of the optical identifications them- 
selves are negligible in comparison. We use the positional 
uncertainty ellipse of FSC92 to find matching counterparts 
of the IRAS sources, instead of using a circular matching 
tolerance (to be discussed in Section 12. 3p . If a galaxy in 
a redshift survey lies within 3 a positional uncertainty el- 
lipse of the IRAS source, we regard it a match. As a re- 
sult of this matching, we find that 8382, 2091 and 10197 
IRAS sources have optical counterparts in the SDSSDR4, 
2dFGRSDRF and 6dFGSDR2, respectively. Some IRAS 
sources (615 for SDSSDR4, 201 for 2dFGRSDRF, and 91 
for 6dFGSDR2) have more than two galaxies within 3 a 
error ellipse. In these cases, we computed the "Likelihood 
Ratio (LR)" (|Sutherland fc Saunderslll992l . to be discussed 
in I2.3|) of each association, and selected a more appropriate 
galaxy with the larger value of LR as an optical counterpart. 
In total, we compiled a list of 19,380 sources by collecting all 
matched sources from our disparate galaxy redshift surveys. 



2.2 Selection Criteria 

We used the following steps in order to identify bona-fide 
ULIRGs and HLIRGs from the list of all matched sources: 

(i) For 19,380 sources in the list of all matched sources, we 
calculated the IR luminosities using the IRAS 60 u.m fluxes 
and the redshifts of galaxies on the assumption that the 
ULIRGs and HLIRGs have M82 starburst Spectral Energy 
Distributions (SEDs). A substantial fraction of the IRAS 
sources have upper limit fluxes at 12 (i.m, 25 [im and 100 
(i.m, so we use this single-band selection to obtain a less 
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Table 1. ULIRGs found in SDSSDR4, 2dFGRSDRF and 6dFGSDR2 



FSC 
NAME 


SDSS 






6dFGS 






2dFGRS 






Final 


ID 


z 


zconf 


ID 


z 


Q 


ID 


z 


Q 


F00050-3259* 














TGS434Z061 


0.285 


4 


3 


F00091-0738 








g0011433-072207 


0.118 


i 








2 


F00091-3905* 








g0011432-384913 


0.254 


3 








2 


F00095-5948* 








g0011586-593133 


0.235 


4 








2 


F00184-3331* 








g0020557-331428 


0.239 


3 


TGS497Z156 


0.239 


4 


2 


F00256-0208 








g0028143-015146 


0.277 


3 








2 


F00285-3140 








g0030593-312445 


0.218 


2 


TGS439Z175 


0.217 


4 


3 


F00318-3137* 








g0034133-312118 


0.285 


4 








2 


F00335-2732 














TGS206Z015 


0.069 


5 


3 


F00406-3127 








g0043032-311050 


0.343 


1 








2 



Column descriptions: (1) The IRAS object name in the FSC92. Asterisks represent the sources that are identified as ULIRGs for the 
first time in this study; (2-4) The identification, the redshift, and the redshift confidence value in SDSSDR4; (5-7) The identification, 
the redshift, and the redshift quality parameter in 6dFGSDR2; (8-10) The identification, the redshift, and the redshift quality 
parameter in 2dFGRSDRF; (11) Finally accepted optical counterpart (1— SDSSDR4, 2-6dFGSDR2, 3— 2dFGRSDRF). 



Table 2. The IRAS data for the Final Sample of ULIRGs 



FSC 


RA 


Dec 




12imi 




25um 




60 um 




100 Mm 


Flux 


NAME 


(J2000) 


(J2000) 




(Jy) 




(Jy) 




(Jy) 




(Jy) 


Qual 


F00050- 


-3259* 


00 07 34.6 


-32 43 03 


0.066 


(0.018) 


0.144 


(0.042) 


0.222 


(0.045) 


0.758 


(0.220) 


1131 


F00091- 


-0738 


00 11 43.3 


-07 22 05 


0.071 


(0.020) 


0.215 


(0.054) 


2.626 


(0.184) 


2.521 


(0.202) 


1232 


F00091- 


-3905* 


00 11 42.3 


-38 49 15 


0.106 


(0.030) 


0.125 


(0.038) 


0.316 


(0.047) 


0.756 


(0.227) 


1131 


F00095- 


-5948* 


00 11 58.8 


-59 31 28 


0.074 


(0.014) 


0.052 


(0.010) 


0.313 


(0.034) 


0.681 


(0.123) 


1132 


F00184- 


-3331* 


00 20 57.7 


-33 14 28 


0.105 


(0.028) 


0.120 


(0.035) 


0.334 


(0.047) 


0.613 


(0.141) 


1132 


F00256- 


-0208 


00 28 14.4 


-01 51 47 


0.108 


(0.031) 


0.323 


(0.090) 


0.602 


(0.060) 


0.611 


(0.147) 


1132 


F00285- 


-3140 


00 31 03.0 


-31 24 18 


0.144 


(0.042) 


0.141 


(0.041) 


0.389 


(0.051) 


0.688 


(0.199) 


1132 


F00318- 


-3137* 


00 34 16.0 


-31 21 04 


0.108 


(0.030) 


0.167 


(0.048) 


0.257 


(0.062) 


0.563 


(0.158) 


1132 


F00335- 


-2732 


00 35 59.2 


-27 15 42 


0.144 


(0.045) 


0.632 


(0.069) 


4.294 


(0.472) 


3.207 


(0.225) 


1332 


F00406- 


-3127 


00 43 03.0 


-31 10 53 


0.060 


(0.016) 


0.091 


(0.025) 


0.717 


(0.057) 


0.994 


(0.169) 


1132 



Column descriptions: (1) The IRAS object name in the FSC92. Asterisks represent the sources that are identified as ULIRGs for the 
first time in this study; (2-3) Right ascension and declination; (4-7) The IRAS flux density (and its uncertainty) at 12(jm, 25um, 
60um and lOOum; (8) The IRAS flux density qualities at each band (1— upper limit, 2— moderate quality, 3— high quality). 



heterogeneous final catalogue (to be discussed in Section 

(ii) We selected the sources for which their 60 u.m fl ux 
qualities are "high" or "moderate" l|Moshir et al.lll992l ) in 
order to restrict our analysis to those sources with reliable 
60 u.m detection^], obtaining 19,335 sources. 

(iii) Of the 19,335 reliable sources, we classified 483 
sources as ULIRG candidates (IR luminosities greater than 
10 12 Lq), and classified 14 sources as HLIRG candidates (IR 
luminosities are greater than 10 13 L©). 

(iv) We then assessed the quality of the redshift for the 
ULIRG and HLIRG candidates. We required a redshift con- 
fidence of ^ 0.65 for the SDSSDR4, a redshift quality of ^ 
3 for the 2dFGRSDRF, and a redshift quality is equal to 3 
or 4 for the 6dFGSDR2. See each data release paper for the 
descriptions of these redshift quality parameters. On this ba- 



4 The flux measurements of high quality, moderate quality or 
upper limits are represented as number of 3, 2, 1 respectively in 
FSC92. 



sis, we secured 333 ULIRG and no HLIRG candidates with 
acceptable redshift qualities. 

(v) We used the NASA/IPAC Extragalactic Database 
(NED) Near Position Search in order to reject chance co- 
incidences of IRAS sources for which optical counterparts 
are already listed in the literature. We rejected 13 ULIRG 
candidates (two sources moved to a list of HLIRGs), and 
added 4 sources (to be discussed in Section [2.3[) . Our red- 
shifts for these four are confirmed in the literature, although 
their redshift qualities do not satisfy our criteria. Finally we 
found 324 ULIRGs and 2 HLIRGs (the final two HLIRGs 
will be discussed in Section T4.3[l , 

Table[T]lists the final sample of 324 ULIRGs with the names 
of the IRAS sources, the identifications of their optical coun- 
terparts, redshifts and redshift quality parameters in each 
galaxy redshift survey. In some cases, more than two dis- 
parate galaxies from different galaxy redshift surveys are 
located within an error ellipse of one IRAS source. There- 
fore, we present the flags indicating which galaxy has been 
adopted as the optical counterpart in the final column of 
Table [1] We present the IRAS FSC92 data in Table [2] 
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Table 3. The Final Sample of ULIRGs 



FSC 


Final 




log(LR) 


Class 


NAME 


z 


(£©) 




L/S 


F00050- 


-3259* 


0.2855 


12.12 


0.5042 


7/V 


F00091 


-0738 


0.1178 


12.30 


0.6345 


7/IV 


F00091- 


-3905* 


0.2535 


12.14 


0.5944 


7/V 


F00095- 


-5948* 


0.2349 


12.06 


0.6102 


7/IV 


F00184- 


-3331* 


0.2387 


12.10 


-0.1065 


7/V 


F00256 


-0208 


0.2770 


12.52 


0.5735 


5/IIIa 


F00285 


-3140 


0.2166 


12.07 


-0.9137 


4/1 


F00318- 


-3137* 


0.2846 


12.18 


0.3252 


7/V 


F00335 


-2732 


0.0686 


12.01 


0.3877 


7/IV 


F00406 


-3127 


0.3425 


12.82 


0.5803 


7/V 



19991: Lawrence et al.l 19991; RIeodouIou et al.l 


1999; 


Stanford et al. 20001; ISanders et al. 


20031; Gotol 


2005; 


Cao et al.ll2006f). We then aDDlied the 


same procedure given 



Column descriptions: (1) The IRAS object name in the FSC92. 
Asterisks represent the sources that are identified as ULIRGs for 
the first time in this study; (2) Finally accepted redshift from 
galaxy redshift survey; (3) Infrared luminosity calculated in 
this study; (4) Likelihood ratio; (5) Morphological class using 
Lawrence classification (L) and Surace classification (S). See the 
description in the end of this section. 



Adopted redshifts, IR luminosities, log (LR) and morpho- 
logical classes (to be discussed at the end of this section) 
are listed in Table [3] 

Figure [1] shows a sample of 3' x 3' finding charts (left 
hand panels) for the final sample of ULIRGs extracted 
from the SDSS gri composite images or the SuperCOS- 
MOS scans of tf survey plates. The contour plot for each 
ULIRG using the SDSS r-band images or the SuperCOS- 
MOS tf images is also presented in Figure [T] (right hand 
panels). Using Figure [T] we classify the interaction types 
of the ULIRGs in the final sample. The ULIRGs are clas- 
sified according to the m odified version o f the interaction 
classification schemes of lLawrence et al. I (Il989l . hereafter 
La wrence classifi cation) (see, e.g.. iFarrah et al] 120011 ) and 
of Suracd |l998|. hereafter Su race classification) (see, e.g., 
IVeilleux. Kim fc Sandersll2002r i , as are summarised in Table 
21 The results of the classification are given in the last col- 
umn of Table [3] The complete versions of Table [T] Table [2j 
Tableland Figure [T] are available only in the electronic is- 
sue. A brief discussion of some ULIRGs is given in Appendix 
A that is available in the electronic issue. 



2.3 Reliability of Identifications 

A variety of definitions of ULIRGs have been used in the 
literature, and depend on the minimum luminosity, the 
adopted spectral energy distributions and the cosmology 
(L 8 _iooo ^ 1 12 L B with Hp =75 km s^Mpc" 1 and 
go = 0.0 for ISanders fc Mirabel Il99(l iKim fc Sandersl 



1 19981 . and IGotol 120051 : 

km s _1 Mpc _1 and go 
L 60 > 10 12 L W 



Leo ^ iir — l,q witn ti p- 

= 0.5 for IClements et~ai1 Il996al ; 
with Hp= 50 km s^Mpc" 1 and fi=l for 
lLawrence et alj Il999f) . In order to check the reliability 
of our selection criteria of ULIRGs, we compiled 636 
known ULIRGs taken mostly from wide area survey 



in Section 12.21 to select ULIRGs in this sample of known 
ULIRGs. Of these known ULIRGs from the literature, 



data ( 


Leech et all 19941; Clowes et all 1995; Murohv et all 


1996: 


Clements et al.l 


1996al; iDuc. Mirabel & Maza 1 19971; 


Kim & Sandersl 1998 


; IClements, Saunders & McMahon 



we found 291 reliable sources in step (ii) of Section 12.21 
A total of 345 known ULIRGs were not found since they 
were neither observed in the redshift surveys of this study 
nor included in FSC92. After the final identification step 
in Section 12.21 we found 131 known ULIRGs in our final 
sample of ULIRGs, and failed to find a total of 160 ULIRGs. 
The reasons our having missed these 160 known ULIRGs 
are categorised as follows; 

(i) Different redshifts (F00415-0737, F01031-2255, 
F01082-2452, F03014-2026, F03193-2224, F03448-2628, 
F08509-1504, F10122+4943, F20087-0308, F20109-3003, 
F20175-4756, F23515-2917): For the sources in this cate- 
gory, the redshifts used in this study are different to (and 
almost always much lower than) those in the literature. 
Therefore, the IR luminosities calculated in this study are 
not large enough to be selected as ULIRGs in this study. 
Some sources have high redshift qualities (F00415— 0737, 
F08509-1504, F10122+4943, F20109-3003), while the red- 
shift qualities for the others are low. The redshift of F01031- 
2255 used in this study is larger than that in the literature, 
but it is not included in the final sample of ULIRGs due to 
its low redshift quality. 

(ii) Different optical counterpart galaxies 
(F03202-0001, F08007+3928, F09346+3911, F12527-0306, 
F14390+6209, F14475+0155, F14546+0338, F15182+3023, 
F21368+1006, F22011+0017, F23051-0100): The optical 
counterpart galaxies of these IRAS sources in this study 
are different from the galaxies in the literature, as are their 
redshifts. 

(iii) Different 60^m fluxes (F00444-1803, 
F09320+6134, F11087+5351, F14170+4545, F14351-1954, 
F14575+3256, F18520-5048): The IR luminosities of these 
sources are not greater than 10 12 L© since the 60 \xm fluxes 
of these sources used in this study (from FSC92) are lower 
than the fluxes used in the literature (mostly from PSC). 

(iv) Marginal ULIRGs (F12495-3414, F13156+0435, 
F14378-2604, F22509-0040): Since the L 60 luminosities 
of all four sources are slightly less than 1O 12 L0, they 
were classified as fu rther marginal/possible ULIRGs in 
lLawrence et al.l (|l999l ). Similarly, their IR luminosities cal- 
culated in this study are in the range of 10 11 ' 84 — 10 1195 Lq, 
therefore they are not classified as ULIRGs in this study. 

(v) Low qualities of redshifts (F00406-3127, 
F10026-0022, F14207-2002, F14485-2434, F22123-2025): 
These sources satisfy the selection criteria up to step (iii) of 
Section \2. 21 However, they fail to satisfy the criteria in step 
(iv) of Section 12.21 due to their low redshift qualities. Al- 
though the redshift qualities of F00406-3127, F14485-2434 
and F22123— 2025 are low, these redshifts are same as those 
quoted in the literature. Therefore we include these three 
sources in the final sample of ULIRGs. F14207-2002 is 
not a know n ULIRG, but its redshift in this study is same 
as that of lAllen et all l|l99ll ). Therefore we also include 
F14207-2002 in the final sample of ULIR Gs. H o wever , 
F10026-0022 was identified as a ULIRG in IGotol (|2005l ) 
using the same SDSS data with low quality of redshift as in 
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Figure 1. The 3'x3' finding charts (left in each panel) and contour plots (right in each panel) of 324 ULIRGs identified in this study. 
The finding charts centred on IRAS source positions are extracted from the SDSS gri composite images or the SupcrCOSMOS scans of 
rp survey plates. North is up, and the East is to the left. The ellipse represents 3<r IRAS positional uncertainty. The optical counterpart 
in the redshift surveys is marked by square indicating the size, orientation and location of the contour plot in the finding chart. The 
IRAS object names are presented in the bottom left corner, and asterisks represent the newly identified ULIRGs in this study. The 
contour plots centred on optical counterparts represent the intensities of ULIRGs in the SDSS r-band images or the SuperCOSMOS r^ 
images. The size of each contour plot is 40 kpc X 40 kpc. The contour levels increase from the sky level (/ s k y ) to the peak value (/ ma x) 
by factors lo' Imax_ J sk y )/8_ ^phe orientation is indicated by the arrows and the thick, horizontal bar represents 5 arcsec in each contour 
plot. 
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Table 4. The interaction classification scheme of ULIRGs used in this study 



Scheme 


Class 


Description 


Lawrence 





Isolated source with no signs of interaction or merging 


classification 


1 


Source with a faint companion (2~4 mag fainter than the source) in the range of 40~200 kpc 




2 


Source with a bright companion (less than 2 mag fainter than the source) in the range of 40^200 kpc 




3 


Source with a faint companion less than 40 kpc 




1 


Source with a bright companion less than 40 kpc 




5 


Source interacting with a companion and showing signs of interaction 




6 


Merger/More than two nuclei in common envelope 




7 


Merger/Single nuclei in common envelope 


Surace 


I 


First approach — Unperturbed and separated galaxies with no signs of interaction or merging 


classification 


11 


First contact - Overlapped disks but no strong bars and tidal tails 




Ilia 


Pre-merger a — Two recognisable nuclei with strong signs of interaction and separated more than 10 kpc 




Illb 


Pre-merger b — Two recognisable nuclei with strong signs of interaction and separated less than 10 kpc 




IV 


Merger - Only one nucleus seen with strong tidal features 




V 


Old merger — No direct signs of tidal features but disturbed central morphology 



this study. Therefore we reject it from the final catalogue of 
ULIRGs. 

(vi) Underluminous ULIRGs (F00090-0054, 
F08112+3039, F08322+3609, F09045+3943, F16134+2919, 
F 21341-0033): These sources are the FIRST-FSC sample 
of I Stanford et all (|2000l ) . They did not use a strict minimum 
infrared luminosity to identify ULIRGs. T he far-infrared 
(FIR) luminosities for these sources in IStanford et al.l 
|2000l ) are in the range of 10 11 ' 27 - 10 11 85 L Q which is 
similar to the minimum FIR luminosity of ULIRGs in 
ISanders fc Mirabell (|l996h . Similarly, these sources are not 
selected as ULIRGs in this study due to our low estimates 
for their IR luminosities (10 11 ' 44 - 10 11 " L ). 

(vii) Different SEDs (F 21219-1757): Th i s sour ce was 
classified as a ULIRG in iKim fc Sanders! (|l99St ) with 
log(Li r /L0)=12.O6. However, the IR luminosity calculated 
in this study is log(L; r /L0)=11.86 which does not satisfy the 
ULIRGs criteria, although the redshift and IRAS flux densi- 
ties us ed in this study are same as those in lKim fc Sanders] 
(1998). This is because we used M82 starburst SED to cal- 
culate the IR luminosity which may not be appropriate for 
low luminosity QSOs such as this. If we use AGN dominated 
SEE0 of NGC 1068 to calculate the IR luminosity for this 
source, we obtain the IR luminosity of log(Li r /L0)=12.14, 
satisfying the ULIRGs criteria. 

(viii) Using upper limit flux (11 5 obj e cts): These 
sources were classified as ULIRGs in iGotol (|2005l) and 
ICao et all (120061 ). Although most IRAS galaxies have up- 
per limit fluxes at 12 u.m and 25u .m (to be discussed in the 
end of this section), IGotol (|2005l ) calculated the IR lumi- 
nosities of IRAS sources by treating the flux upper limits 
as detections. Similarly, three sources for which 100|J.m flux 
densities are "u pper limits" were found to be ULIRGs in 
ICao et all (|2006l ) by assuming the 100 (im upper limit is ob- 
tained. These calculations can clearly result in overestimates 
for the IR luminosities of IRAS galaxies. These sources are 



5 If we use AGN dominated SED of NGC 1068 to calculate the 
IR luminosities for ' Warm' ULIRGs (/2s //60 > 0.2) in the final 
sample, the IR luminosities for those sources will be increased by 
82 per cent. 




Figure 2. The 2'x2' finding chart of F12265+0219 extracted 
from the SDSS gri composite image. The finding chart centred 
on IRAS source position. North is up and East is to the left. 
The ellipse represents 3a positional uncert ainty, a nd the square 
denotes the optical counterpart found in I Gotol ll2005t) . SDSS 
J122907.3+020246. The star like object within the ellipse is 3C 
273. 



not selected in our study as ULIRGs due to their low IR 
luminosities of 10 10 ' 60 - 10 11 ' 99 L . 

In summary, all 160 known ULIRGs that are not iden- 
tified in this study have appropriate reasons not to be se- 
lected as ULIRGs except for only one source (F21219— 1757). 
Therefore we conclude that our ULIRGs selection criteria 
are reliable and good enough to identify nearly all known 
ULI RGs 

IGotol (|2005l ) performed a cross-correlation of FSC92 
with the spectroscopic catalogues of galaxies in the SDSS 
Third Data Release (SDSSDR3), and identified 178 ULIRGs 
and 3 HLIRGs. Since the SDSSDR4, which we used for 
our cross-correlation in this study, contains galaxies of 
SDSSDR3, we can compare our iden tificati o n pro cedure of 
ULIRGs and HLIRGs with that of iGotd (|2005l ). Firstly, 
we accept the correlation if a galaxy in SDSSDR4 lies 
withi n 3 a position al uncertainty ellipse of the IRAS source, 
while IGotol l|2005l ) treats a neighbouring optical galaxy as 
a match if the galaxy in SDSSDR3 lies within 20 arc- 
sec of the IRAS source without regard to the IRAS er- 
ror ellipse. This can cause many problems. For example, 
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SPSS J122907.3+020246 was identified as a HLIRG in lGotol 
|2005l ). In Figure [2] we present the SDSS gri composite im- 
age of F12265+0219, which was regar ded as the I RAS coun- 
terpart of SDSS J122907.3+020246 in lGotol (|2005T ). It is seen 
that SDSS J122907.3+020246 lies outside the 3a error el- 
lipse and a bright star-like object lies within the ellipse. The 
star-like object within the ellipse is the famous quasar 3C 
273, which is known to be associated with F12265+0219. 
Since 3C 273 wa s not i nclude d in the spectroscopic targets 
of the SDSS and lGotol (|2005T l used circular matching toler- 
ance of 20 arcsec, the mismatch between F12265+0219 and 
SDSS J 122907.3+020246 was included in his sample. Our 
cross-correlation based on the positional error ellipse, and 
the NED Near Position Search in step |(v)| of Section ^. 2l help 
to avoid this kind of mismatching. Secondly, there is a differ- 
ence i n the calculation of infrared luminosity between iGotol 
(2005) and this study. The forme r used the same function 
as used by iKim fc Sander s (1998) for the calculation of IR 
luminosity using IRAS flux directly from FSC92. However, 
a significant fraction (81 per cent of "high" or "moderate" 
60 u.m sources for /12 fluxed and 74 per cent for /2s fluxes) 
of the IRAS sources only have upper limits for /12 and /2s 
fluxes. Therefore, assuming upper limits are obtained for the 
measurements at /12 and /2s in calculating the IR luminos- 
ity can lead to overestimates of the IR luminosity, although 
the uncertainty in the IR luminosity is only a few percent 
l|Kim fc Sanderslll998h . 

2.4 Reliability and Completeness of the final 
sample 

A final sample of ULIRGs selected by cross-correlation based 
on the position alone may contain spurious sources due to 
the chance presence of a galaxy within the IRAS 3 a ellipse. 
In order to estimate the probability of the "true" associa- 
tion between IRAS source s and the optical counterparts , we 
compute likelihood ratio (|Sutherland fc SaundersHl992T l for 
each association. The likelihood ratio defined by the ratio 
of the probability of a true association to that of a chance 
association ifl 



LR : 



Q« m)exp(-r 2 /2) 
27ro"i<T2n(^ m) 



(1) 



where m) is the local surface density of objects brighter 
than the candidate. The "normalised distance" r is given by 



(2) 



where the di, di are the positional differences along the two 
axes of an error ellipse for an IRAS source, and ci, 02 are the 
lengths of these axes. Since the positional uncertainties of 
galaxies in SDSS, 2dFGRS, and 6dFGS are negligible com- 
pared with those of IRAS sources, we define a as the length 
of the error axes of the IRAS sources. m) is a multi- 
plicative factor which is the a priori probability that a "true" 
optical counterpart brighter than the magnitude limit exists 
in the association, and for simplicity we set Q = 1 in this 
study. 




Figure 3. Reliability of our cross-correlation for SDSSDR4 (dot- 
ted line), 2dFGRSF (dashed line) and 6dFGSDR2 (solid line) as a 
function of LR. The vertical lines represent the critical LR values 
of reliable identification for each redshift survey data - above the 
critical LR values, the reliabilities are 80 %. 



We compute n(^ m) using photometric sample of galax- 
ies within 3 a error ellipse, 



n(^ m) 



Nj£ m) 

97T<T1(T2 



(3) 



where m) represents the number of galaxies of which 

magnitude^] are less than or equal to that of a candidate. 
We then obtain the likelihood ratio for our sample 



LR : 



9exp(-r 2 /2) 
2iV(^ m) ' 



(4) 



We compute LR values for all 19,380 sources in the list of all 
matched sources of Section [2.11 using the photometric sam- 
ple of galaxies in each redshift survey. In order to calculate 
the reliability of association using the LR values, we perform 
random associations by offsetting the IRAS source positions 
by ~ 30', and recompute LR values for e ach random associa - 
tion following lLonsdale et ail (|l998r ) and lMasci et all (|20Q]J) . 
Using the distribution of LR values for true associations and 
random associations, the reliability of each association with 
a given LR is defined by 



R(LR) = 1 



N„ 



■ (LR) 



JVtrue(LR) 



(•») 



where iVt ruc (LR) and iV ran dom(LR) are the number of true 
and random associations with a given LR. In Figure [3] we 
present the reliability of our cross-correlation for each red- 
shift survey as a function of LR. Above log (LR) ~ —0.44 
(SDSSDR4), -0.39 (2dFGRSDRF), -1.29 (6dfGSDR2), the 
reliabilities are 80%. We show the LR distribution of our 



6 The quantities /12, /25, /60, and /100 represent the IRAS flux 
densities in Jy at 12^im, 25p_m, 60^m, and 100|^m, respectively. 

7 We assume that a positional error of IRAS source is a Gaussian. 



8 We use SDSS Petrosian r-band magnitude for the galaxies in 
SDSSDR4, and bj band magnitude for the galaxies in 2dFGRS- 
DRF and 6dFGSDR2. 
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Figure 4. Distribution of LR values for the final sample of 324 
ULIRGs found in SDSSDR4 (shaded region), 2dFGRSDRF (filled 
region), and 6dFGSDR2 (solid line). The vertical lines represent 
the critical LR values of reliable identification (the reliability of 
~ 80 %) as shown in Figure [3] 



ULIRG sample in Figure [4] along with critical LR values of 
reliable identification for each redshift survey. If we apply 
the critical LR values of 80% reliability, then 24 out of 126 
(SDSSDR4), 2 out of 27 (2dFGRSDRF), and 1 out of 171 
(6dFGSDR2) ULIRGs can be rejected from a final sample 
of ULIRGs. 

The reliability of our LR-based identifications is con- 
firmed by the NED Near Position Search in Section 12.21 Of 
the 11 objects that were rejected from our final sample of 
ULIRGs because the optical association differs from that in 
NED, eight have lower LR values than the critical values. 
Although three objects have larger LR values than the criti- 
cal values, there are more appropriate galaxies in NED (the 
larger value of LR) that would result in the galaxies not 
being ULIRGs. 

The spectroscopic target selection in a redshift survey 
introduces incompleteness into our final sample of ULIRGs. 
Since the target selection functions differ in our three red- 
shift surveys (a short discussion is presented in Section T4.ip 
and are beyond the scope of this paper, we only estimate 
the incompleteness originated from this study. Therefore, 
we measure the completeness relative to all FSC92 ULIRGs 
that are identified in these surveys, and not for example rel- 
ative to a volume-limited ULIRG survey or to FSC92 as a 
whole. 

There are three factors which make our sample incom- 
plete. Firstly, since we used 3 a error ellipse for our cross- 
correlation, some ULIRGs outside 3 a error ellipse might be 
missed. If we assume an error ellipse to be a Gaussian distri- 
bution, the incompleteness originated from this term would 
be about 1 per cent. Secondly, some ULIRGs for which red- 



part galaxy are within the error ellipse of one IRAS source, 
the wrong galaxy may be selected as the optical counterpart 
of the IRAS source. In order to quantify the incomplete- 
ness introduced by those factors, we compute the ratio of 
ULIRGs that are surveyed in SDSSDR4, 2dFGRSDRF and 
6dFGSDR2 to those that are missed from our final sam- 
ple of ULIRGs using previously known ULIRGs. Of the 636 
known ULIRGs in Section 12.31 we consider only the 140 
known ULIRGs that are expected to be found in this study. 
This sample consists of ULIRGs that are located within the 
redshift survey region and are found in FSC92. Addition- 
ally, their IR luminosities calculated using our method are 
greater than 10 12 L©. In the result of identification proce- 
dure in Section 12.11 and 12.21 133 ULIRGs are found in our 
final sample of ULIRGs, and 7 ULIRGs are missed. 

Of the 7 missed ULIRGs, three IRAS sources (F03202- 
0001, F14546+0338 and F22011+0017) have more than two 
galaxies including a known optical counterpart within their 
error ellipses. Our selection procedure is to choose the most 
likely counterpart (i.e. the galaxy with the larger value of 
LR), but this selects different galaxies to the known opti- 
cal counterparts for these three IRAS sources. Therefore, 
the incompleteness introduced by our selection using LR 
values among several counterpart candidates is about 2.1 
per cent (3 out of 140), although all the optical counter- 
parts that we select are brighter or closer to the IRAS posi- 
tion than the known optical counterparts. In addition, one 
ULIRG (F09111-1007) is missed, because the optical coun- 
terpart in 6dFGSDR2 is outside 3 a error ellipse, although 
the association is already confirmed in the literature. Three 
ULIRGs (F00406-3127, F14485-2434, and F22123-2025) 
are missed because of their low redshift qualities as shown 
in Section [2. 31 although the redshifts quoted in the redshift 
surveys are same as those in the literature. Therefore, the in- 
completeness introduced by our cross-correlation algorithm 
using 3 a error ellipse is about 0.7 per cent (1 out of 140), 
and the incompleteness introduced by our selection using 



the redshift quality parameter in the step (iv) of Section 



shift qualities are unreliable (step (iv) in Section ^. 2p are re- 



jected, although their redshifts are not wrong. Thirdly, when 
more than two galaxies including a known optical counter- 



I2.2l is about 2.1 per cent (3 out of 140). In total, the incom- 
pleteness of the final sample of ULIRGs originated from this 
study is about 5 per cent (7 out of 140). 



3 RESULTS 

We have found 324 ULIRGs (126 in SDSSDR4, 27 in 2dF- 
GRSDRF, and 171 in 6dFGSDR2), 190 of which are newly 
discovered. We have therefore increased the number of 
known ULIRGs by about 30 per cent. 

In Figure [5] we present the redshift distribution for the 
final sample of 324 ULIRGs identified in this study, and com- 
pare it to that of the 636 previously known ULIRGs. The 
redshifts in our final sample run from z=0.037 to z=0.517 
and the median value is z=0.223, which is larger than 
the median redshift z=0.184 of the 636 previously known 
ULIRGs. 

In Figure [U we plot the IR luminosities of the ULIRGs 
identified in this study against their redshifts. In the upper 
and right panel, we present the redshift and luminosity dis- 
tribution of ULIRGs found in each redshift survey, respec- 
tively. The ULIRGs identified in SDSSDR4 are found from 
low z to high z and from low luminosity to high luminosity 
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Figure 5. Redshift distribution of the 324 ULIRGs identi- 
fied in this study (shaded region) compared to the 636 known 
ULIRGs (dashed lines), 134 known ULIRGs in the final sample 
in this study (filled region) and known ULIRGs plus all identified 
ULIRGs in this study (solid lines). 
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Figure 6. The IR luminosity of identified ULIRGs in this study 
versus redshift. The ULIRGs from the SDSSDR4, 2dFGRSDRF 
and 6dFGSDR2 are represented by squares, triangles, and pluses, 
respectively. In the upper and right panel, the redshift and lumi- 
nosity distribution of ULIRGs found in each redshift survey are 
shown, respectively. 

(~ 1O 12 ' 85 L ). The ULIRGs in 6dFGSDR2 show similar dis- 
tribution with those in SDSSDR4, although median redshift 
of galaxies in 6dFGSDR2 is smaller than that of galaxies in 
SDSSDR4. In contrast, in spite of relatively larger median 
redshift of galaxies in 2dFGRSDRF compared to that of 
6dFGSDR2, no ULIRGs in 2dFGRSDRF are found beyond 
z~0.35. This might be due to the different selection crite- 
ria for spectroscopic targets (see each data release paper for 
more detail). Figure [6] shows that the upper envelope (shown 
by a dotted line) increases from log (Li r /L@) = 12.45 to 
12.95 as redshift increases. The lower envelope is due to the 
survey detection limit. 




12.0 12.2 12.4 12.6 12.8 13.0 

iog(iVU> 

Figure 7. The luminosity distribution for the final sample of 
324 ULIRGs is plotted in the top panel. Below this panel, the 
luminosity distributions of ULIRGs in the redshift ranges 0.00 ^ 
z < 0.15, 0.15 < z < 0.2, 0.2 < z < 0.25, 0.25 < z < 0.35, and 
0.35 ^ z are shown. 



The luminosity distribution of the final sample is plot- 
ted in the Figure [7] For the final sample of ULIRGs (top 
panel), there are fewer ULIRGs at high luminosities than at 
low luminosities. The median luminosity value is 10 12 ' 17 Lq. 
Two interesting features are seen in the luminosity variation 
with redshift in Figure [7J Firstly, the number ratio of low 
luminous ULIRGs (L ir < 1O 12 2 L ) to all ULIRGs in red- 
shift bin increases with redshift up to z=0.25 (due to the 
detection limit seen in Figure there are few low luminous 
ULIRGs beyond z=0.25). Secondly, no ULIRGs that have 
IR luminosities of log (Li r /Lo) > 12.60 are found in the 
lower redshift (z<0.25). This is almost certainly due to the 
evolving luminosity function. A calculation of the evolving 
luminosity function is beyond the scope of this paper. 



4 DISCUSSION 

4.1 Comparison of ULIRG subsamples 

Since we use three different redshift survey data to search 
for ULIRGs, different target selection functions in each red- 
shift survey make our final sample of ULIRGs inhomoge- 
neous. The majority of spectroscopic targets in each red- 
shift survey are selected from magnitude limited samples^], 

9 The primary sample (~ 70%) of 6dFGS are galaxies with 
Ktot <12.75 and one additional sample is optically selected galax- 
ies with bj < 16.75 (~ 5%). Most spectroscopic targets of 2dF- 
GRS are galaxies with bj < 19.45, and the main galaxy sample 
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Figure 8. Absolute bj magnitude of ULIRG subsamples against 
redshift. The ULIRGs from the SDSSDR4, 2dFGRSDRF and 
6dFGSDR2 arc represented by squares, triangles, and plusses, re- 
spectively. The dotted and dashed line indicate the bj magnitude 
limit of 6dFGSDR2 and 2dFGRSDRF, respectively. 
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Figure 9. Infrared colour-colour diagram. Log(/2s//6o) i s plot- 
ted against log(/6o//ioo) for the ULIRG subsamples. The 
ULIRGs from the SDSSDR4, 2dFGRSDRF and 6dFGSDR2 are 
represented by squares, triangles, and plusses, respectively. 



but some additional spectroscopic targets are included de- 
pending on the redshift survey. In order to see the effects 
of the magnitude limit of each survey on our ULIRG sub- 
samples, we plot the absolute bj magnitude against red- 
shift of ULIRG subsamples in Figure [5] A transformation 
from SDSS photo metry to bj m a gnitud e, and k-correction 
is done following iNorberg et al.l (|2002l ). Since the major- 
ity of the spectroscopic targets in 6dFGS are K-band se- 
lected, most ULIRGs found in 6dFGSDR2 are fainter than 
the bj magnitude limit (bj = 16.75) shown by the dot- 
ted line. The majority of ULIRGs found in 2dFGRSDRF 
are above the magnitude limit (bj = 19.45). The bj mag- 
nitude limit for the spectroscopic sample of SDSSDR4 is 
not clearly defined since the main spectroscopic sample of 
galaxies are Petrosian r-band selected sample. However, bj 
magnitude limit for the main galaxy sample of SDSSDR4 is 
similar to that of 2dFGRSDRF, and most galaxies fainter 
tha n the bj magnitude limit are "L uminous Red Galaxies" 
( see lAdelman- McCarthy et al.ll2005l for more detail). It ap- 
pears that a significant fraction (37%) of ULIRGs found in 
SDSSDR4 are fainter than the magnitude limit of 2dFGRS, 
implying that they are Luminous Red Galaxies. 

In order to compare the infrared properties of our 
ULIRG subsamples associated with each redshift survey, we 
plot the infrared colour of log(/2s//6o) against log(/6o//ioo) 
(a detailed discussion on the infrared colours of ULIRGs will 
be given in Section [4. 2p . Although the majority of ULIRGs 
in the range of log(/6o//ioo)< —0.4 or in the range of 
log(/2s//6o)> —0.7 are IRAS sources with flux quality of 
1 at lOOum or at 25um (see Figure [TT1 and [T2 ^l . there are no 



of SDSS are galaxies with Petrosian r-band magnitudes brighter 
than r=17.77. 



striking differences in infrared colours among the ULIRG 
subsamples associated with each redshift survey. 

In conclusion, the sample of ULIRGs found in the three 
redshift surveys are not a simple magnitude limited sample 
due to the diverse target selection functions of each redshift 
survey. 



4.2 Infrared Colours 

Infrared colours such as log(/i2//eo), log(/2s//6o) and 
log(/eo//ioo) have long been used for classifying infrared 
sources. Infrared-bright stars have SEDs peaking typically 
around 12 urn, and can be distinguishe d from galaxies us - 
ing the colour of log(/i2//6o) (see, e.g.. lCohen et alJll987t ). 
In addition, log(/2s//6o) has been used to classify ULIRGs 
into ' Warm' and 1 Cool' systems. It has been suggested 
that ULIRG may evolve from ' Cool' system to ' Warm' sys - 
tem dSanders et al.lll988bl; IVeilleux. Kim fc Sanders! 120021 ). 



ISoifer fc Neugebauerl 1 1991 ) showed that for luminous in- 
frared galaxies including ULIRGs in the IRAS Bright 
Galaxy Survey, the mean log(/eo//ioo) colour increases with 
increasing IR luminosity. 

We present the IR luminosity against infrared colour 
of log(/i 2 // 6 o) in Figure [TO] against log(/2 5 //6o) in Fig- 
ure [11] and against log(/6o//ioo) in Figure [121 for the fi- 
nal sample of 324 ULIRGs. The infrared- bright stars are 
know n to have log(/i2//6o)> (see, e.g., iKim fc Sanders! 
1998). All ULIRGs in the final sample have infrared colours 
of log(/i 2 // 6 o)< 0. 

In Figure [TT] we plotted the boundary lines for the clas- 
sification of normal galaxies, S e yfert 2 galaxies, and QSOs 
suggested by iNeff fc Hutchingsl ( 19921 ) . If we restrict our at- 
tention to the ULIRGs with flux qualities of 2 or 3 (filled 
circles), the number ratio of these objects is N(Galaxy) : 
N(Sy2) : N(QSO) = 32 : 7 : 29. If we use the classification 
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Figure 10. The infrared luminosity against the infrared colour 
l°g(/i2//6o) f° r tti e final sample of 324 ULIRGs in this study. 
The ULIRGs with flux upper limits at 12 (am are represented by 
open circles with arrows indicating the sense of the limit, and the 
ULIRGs with high or moderate flux qualities are represented by 
filled circles. The dotted line represents the boundary between 
galaxies and stars. 
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Figure 11. As Figure 1101 except using the infrared colour 
l°g(/25//6o)- The open circles and arrows represent the ULIRGs 
which have upper limits at 25 um, and the filled circles denote 
the ULIRGs which have high or moderate flux qualities at 25 |^m. 
The dotted lines represent the boundaries among normal galaxies, 
Seyfert 2 galaxies, and QSOs. 



Figure 12. As Figure 1101 except using the infrared colour 
l°g(/6o//ioo)- The ULIRGs with flux upper limits at 100 um 
are represented by open circles with arrows indicating the sense 
of the limit, and the ULIRGs with high or moderate flux qualities 
are represented as filled circles. The open squares are at the mean 
colour and the central luminosity of ULIRGs in the three luminos- 
ity bins (12.0 ^ log(L ir /L e ) < 12.3, 12.3 ^ log(L ir /L Q ) < 12.6, 
and 12.6 SC log(Li r /L,Q) < 12.9). The vertical errorbars de- 
fine the limiting luminosities of the bins, and the horizontal 
errorbars represent standard deviation of log(/eo//loo) colours 
in each bin. The dotted line represents the s election criteria 
(los( /6o//ioo)> -0.3) of 1 Jy sample of ULIRGs llKim fc Sanders! 
11993) . 



scheme of ' Warm' ULIRGs (iV/eo > 0.2, or log(/ 25 //6o)> 
-0.7) and 'Cool' ULIRGs (iW/eo < 0-2 , or log(/ 25 // 6 o)< 
—0.7) suggested bv lSanders et alJ l)l988bl ). the number ratio 
for the sample becomes N(Cool) : N(Warm) = 38 : 30. This 
ratio is much smaller t han that of 1 Jy sam ple (N(Cool) 
: N(Warm) = 90 : 25, iKim fc Sanders! Il998l ). This differ- 
ence might be related to the different selection criteria for 
ULIRGs, but is likely to be strongly affected by the exclusion 
of ULIRGs which have upper limit flux at 25um. 

In Figure \TZ\ the log(/eo//ioo) colours of ULIRGs with 
flux quality of 2 or 3, range from —0.80 to +0.22 and 
have a mean of —0.19. This mean value i s low compared 
to th ose of other ULIRG samples (see, e.g.. lKim fc Sanders! 
Il998l and references therein). In addition, 21 per cent (42 
out of 203) of ULIRGs with flux quality of 2 or 3, have 
colours of log(/6o//ioo) less than —0.3, while no ULIRGs 
in 1 Jy sample have colours of log(/eo//ioo) less than 
—0.3 due to their selection criteri a (log(/6o//ioo)> —0.3) 
of ULIRGs (|Kim fc Sanders! 1 19981 ). We plotted the mean 
colours in each luminosity bin by open squares to investi- 
gate any dependence of the colour on the IR luminosity. 
It appears that there are no significant IR luminosity de- 
pendence of log(/6o//ioo) in our sample (up to 1O 12 ' 6 L0), 
though our colours are still warmer than those typical in less 
IR-luminous galaxies. 

We compare the IR luminosity and the redshift distri- 
bution of ULIRGs with flux qualities of 1 (open circles in 
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Figure 1101 111! and I12[) to those with flux qualities of 2 or 
3 (filled circles). The mean IR luminosity for ULIRGs with 
flux qualities of 1 at 12 [im is not different from that for 
ULIRGs with flux qualities of 2 or 3 as seen in Figure 1101 
Similar results are found for ULIRGs at 25 urn and at 100 
u.m. The median redshift (z=0.230) for ULIRGs with flux 
qualities of 1 at 12 [im is larger than that (z—0A63) for 
ULIRGs with flux qualities of 2 or 3 at 12 (im. Similarly, 
the median redshifts for ULIRGs with flux qualities of 1 at 
25 \im (z=0.236) and at 100 [im (z=0.262) are larger than 
those for ULIRGs with flux qualities of 2 or 3 at 25 u.m 
(1=0.148) and at 100 ^m (1=0.196). 



4.3 Hyperluminous Infrared Galaxies 

Of the 14 HLIRG candidates satisfying selection criteria 

in Section 12.21 all were ex- 



of HLIRGs up to step (iii 



eluded in futher steps. However, we identified two HLIRGs, 
F01044-4050 and F09105+4108, in the course of NED Near 



Position Search for ULIRGs candidates at step (v) in Sec- 



tion 12.21 We list the 2 HLIRGs identified in this study in 
Table [5] Column (1) lists the IRAS object name in the 
FSC92. Columns (2) and (3) list the J2000.0 source posi- 
tion. Column (4), (5), (6) and (7) give the IRAS flux den- 
sities (and their errors) at 12u.m, 25u.m, 60u.m and 100u.m, 
respectively. The IRAS flux density quality at each band is 
given in Column (8). Column (9) lists the finally accepted 
redshift and Column (10) gives the IR luminosity calculated 
in this study. Log (LR) is given in the final Column. In Fig- 
ure [T31 we present 3'x3' grayscale images and contour plots 
extracted from the SuperCOSMOS scans of tf survey plates 
(left panel), as well as those from the SDSS gri composite 
image and the SDSS r-band image (right panel) for the fi- 
nal sample of two HLIRGs. The identification procedure and 
properties of the HLIRGs are given in Appendix B of the 
electronic issue. 



5 SUMMARY 

We present a new sample of ULIRGs and HLIRGs found by 
cross-correlating the IRAS sources in FSC92 with the spec- 
troscopic samples of galaxies in the SDSSDR4, 2dFGRSDRF 
and 6dFGSDR2. Our primary results are summarised below: 

(i) We have identified 324 ULIRGs including 190 newly 
discovered ULIRGs in the regions of the sky covered by the 
SDSSDR4, 2dFGRSDRF and 6dFGSDR2. We increase the 
number of catalogued ULIRGs by about 30 per cent. 

(ii) The reliability of the cross-correlation is estimated us- 
ing the likelihood ratio method. We compute the likelihood 
ratio of each association for our sample of ULIRGs. The 
completeness of the final sample has been estimated using 
previously known ULIRGs in the redshift surveyed region. 
The incompleteness introduced by our identification proce- 
dure due to the cross-correlation using 3 a error ellipse, and 
the selection using LR value and using the redshift quality 
parameter is estimated to be about 5 per cent. 

(iii) The redshifts in our final sample run from z=0.037 
to z=0.517 and the median value is 2=0.223, which is larger 
than that (z=0.184) in previous ULIRG samples. 



(iv) Two HLIRGs, F01044-4050 and F09105+4108, are 
found in the course of NED Near Position Search of ULIRGs 
candidates. 
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